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ABSTRACT 

v— In~this— report  _an  approach  is  proposed  for  the  construction  of  Liapunov 
functions  for  certain  types  of  second  and  third-order  nonlinear. systems  .  If 
the  system  is  described  by  a  vector  differential  equation^  x-^F~(->k«..Jk..„;r.  -  k  ll 


w«e-re--k-j--y~"k^’"=,'--k^'are-“param'eters-oT.the-sy.s.tem,  a  Liapunov  function 
V  =  x'  Px-  fP  is  a^eonstant- matrix)  which  ensures  stability  of  the  linear  system 

_  o,- 

for  all  values  of-—k~.<-k^<»kf  is  called  a  common  Liapunov  function  (.  CLF  )  in 
the  given  range.  .  While  it  may  prove  difficult  to  determine  such  a  CLF  ,  a 
Liapunov  function  rtiay  be  selected  to  ensure  the  stability 

of  the  linear  system  over  the  entire  range  of  the  parameters.  Under  certain 
conditions,  this  Liapunov  function  may  be  easily  modified  for  use  as  a  Liapunov 


£  AS 


function  for  a  nonlinear  system  in  which  the^kifa^  are  functions  of  the  state 
variables . 

Using  this  approach,  sufficient  conditions  are  determined  for  the 
stability  of  a  differential  equation  of-th«-farm"-v^-”~ 

,  rj 

in  terms  of  the  bounds  on  certain  functionstede-rived~from  f,  g  and  h  ? 

IV 
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STABILITY  OF  CERTAIN  NONLINEAR  DIFFERENTIAL  EQUATIONS 
USING  THE  SECOND  METHOD  OF  LIAPUNOV 

by 

Roger  M.  Goldwyn  and  Kumpati  S.  Narendra 
Division  of  Engineering  and  Applied  Physics 
Harvard  University,  Cambridge,  Massachusetts 

I.  INTRODUCTION 

In  recent  years  considerable  interest  has  been  shown  in  the  study  of 
the  stability  of  dynamical  systems  using  the  second  method  of  Liapunov. 
Using  this  approach,  the  stability  of  a  dynamical  system  is  assured  by  the 
determination  of  a  positive  definite  function  V(x)  whose  derivative  with 
respect  to  time  is  negative  definite.  The  conventional  approach  is  one  of  cut 
and  try  where  a  positive  definite  function  V [x)  is  chosen  as  a  candidate  for 

a  Liapunov  function  and  the  parameters  of  the  system  are  adjusted  to  make 

• 

V(x.)  negative  definite.  While  systematic  ways  of  generating  Liapunov 
functions  do  exist  for  linear  time-invariant  systems  (yielding  Routh-Hurwitz 
conditions  implicity),  the  generation  of  such  functions  for  nonlinear  systems 
is  more  difficult.  Recently,  the  "variable  gradient  method"  [Ref.  3  ]  was 
introduced  as  a  systematic  way  of  generating  Liapunov  functions.  While  this 
approach  may,  in  principle,  be  used  for  any  order  system,  in  practice,  even 
for  a  third-order  system,  it  becomes  fairly  involved,  since  no  set  procedure 
exists  for  the  choice  of  the  elements  of  the  gradient  to  satisfy  the  above 
conditions . 


1- 
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In  this  report  a  method  related  to  the  concept  of  a  common  Liapunov 
function  for  a  linear  system  is  presented  for  obtaining  Liapunov  functions  for 
nonlinear  systems  of  low  order.-- Since, for  a  linear  system,  it  is  possible  to 
derive  necessary  and  sufficient  conditions  for  stability,  the  present  approach 
relates  the  nonlinear  problem  to  a  corresponding  linear  problem  for  which 
stability  conditions  are  well  known.  If  a  linear  system  is  described  by  a  set 
of  differential  equations  &=Fx  ,  V(x)  =  x'Px  is  a  Liapunov  function  if 
Q  is  positive  definite  where  V{x)  =  xl[F' P  +  PF]x=-x'Qx  .  While 
every  positive  definite  P  need  not  yield  a  positive  definite  Q  ,  the  solution 
of  the  equation 

-  F'P  +  P  F  =  -  Q 

for  any  positive  definite  Q  must  yield  a  positive  definite  P  if  the  system  is 
stable.  This  important  result  is  used  throughout  the  development  of  the 
present  report. 


II.  COMMON  LIAPUNOV  FUNCTIONS 
Consider  the  linear  system  of  equations 

x  =  Fx  (1) 

where  F  is  a  constant  (n  x  n)  matrix  and  x  an  ( n  x  1  )  vector.  Let  F 

r\ 

be  dependent  linearly  on  a  parameter  k  where  k  :is  assumed  to  be  in  the 
range 

k  <  k  <  k 

If  V(x)  =  x'P_x  where  P  is  a  positive  definite  (or  >  0)  matrix  is  considered 
as  a  candidate  for  a  Liapunov  function 
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where 

-  Q  =  F1  P  +  PF  .  (3) 

If  Eq.  1  is  asymptotically  stable  in  the  given  range,  then  for  any  symmetric 
Q  >  0  the  unique  solution  of  Eq.  3  for  P  yields  a  P  >  0  .  Since  F  =  F(k)  , 
the  resulting  P  will, in  general, be  a  function  of  k— d.  e.  ,  P  =  P(k)  .  Eq.  3 
may  be  considered  as  a  mapping  of  the  set  JQ  l  of  positive  definite  matrices 

II  '  '  * 

into  the  set  <P{k)S  of  positive  definite  matrices.  If  there  exists  a  matrix  P 

1  1  f  | 

which  is  a  member  of  ^P(k)V  such  that 

Q*  (k)  =  -[F'(k).P*  +  P*  F  (k)  ]  (4) 

is  positive  definite  for  all  k  in  the  given  range,  then 

V*  (x)  =  x'  P  x  (5) 

is  a  Liapunov  function  for  Eq.  1  and  will  be  termed  a  common  Liapunov 
function  (  CLF )  in  the  range  (  k  <  kj<  k  )  [Ref.  1]  . 

A.  CLF  for  First-Order  System 
Consider  the  first-order  system  shown  in  Fig.  1  . 


Figure  1.  First-Order  System 
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The  system  may  be  described  by  the  equation 

x  =  ax  -  kx 


Consider  the  function 


(6) 


v  2 
V  =  px 

Then 


(7) 


V  =  2px  [  ax  -  kx]  =  2p  (a  -  k).  x2  '=  -qx2  *  (8) 

It  is  clear  that  Eq.  6  is  asymptotically  stable  for  k  >  a  or,  for  any  q  >  0  , 
p  of  Eq.  7  is  >  0  in  the  region  k  >  k>a  .  In  order  to  obtain  a  CLF  in 

the  given  range,  p  must  be  chosen  independent  of  k  . 

Choosing 


p  =  (k  -  a)  (9) 

V  =  (k  -  a)x  is  a  CLF  for  the  system  in  the  range  k  >  a  since 

V  =  2  (k  -  a )  (a  -  k )  x2  =  -  qxZ  q>0  .  (10) 

If  k  =  k ( x )  ,  V  =  (k  -  a)x2  may  still  be  considered  as  a  CLF  if 
k  =  Min  [k  ( x )  ]  ,  and  k  >  a  since 

V  =  .2  (k  -  a)  (a  -  k(x))x2  <  0  •  (11) 

The  nonlinear  equation  is  asymptotically  stable  in  the  large  for  k  >  a  . 

This  simple  example  indicates  how, for  certain  nonlinearities, a  CLF 
for  a  linear  system  may  also  be  used  as  a  Liapunov  function  for  the  nonlinear 
system. 
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B.  CLF  for  Second-Order  System 
The  simple  example  of  the  preceding  section  was  presented  to  indicate 
the  first  step  in  an  approach  that  might  hopefully  yield  a  logical  procedure 
for  the  derivation  of  results  for  a  third-order  system.  The  difficulties  of  such 
an  approach  become  obvious  even  when  considering  a  simple  second-order 
system  (Fig.  2). 


f(x) 

1 

X 

? 

- ^ - 

2 

s  +  as  +•  b 

— 

- — 

Figure  2.  A  Second-Order  System 


The  second-order  system  shown  in  Fig.  2  may  be  represented  by  the 
set  of  equations 


X1  =  x2 

x2  =  -  k  (  Xj  )  Xj  -  ax 


(12) 


where  k  ( x^  )  =  b  +  f ( ) 

F 


In  terms  of  Eq.  1 

r  o 


-  k ( Xj  )  -  a 


(13) 


When  k  is  treated  as  a  parameter,  the  range  of  interest  is  k  >  0  .  [  For 

k  =  0  one  of  the  eigenvalues  is  zero  so  that  the  system  is  stable,  but  not 
asympotically  stable.] 
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For  positive  definiteness  of  Q  (k)  it  is  necessary  that 
(i)  p12  ,  ap22  -  p12  >  0 

(ii)  det  QJk)  =  q(k)  >  0. 

q  (k)  =  -k2p222  +  4kp12  (ap22  -  p12  )  . 


(18) 

(19) 
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It  should  be  observed  that  it  is  impossible  to  get  a  CLF  of  the  assumed 
quadratic  type  for  all  k  >  0  ,  since  for  some  value  of  k,  q(k)  must  become 
negative.  The  maximum  value  of  k  for  which  a  CLF  may  be  determined  is 
given  by 

r-  4p12<aP22-Pl2> 

k  - - 2 -  ‘  (20) 

p22 


Any  matrix  P  which  satisfies  conditions  Eqs.  16  ar.d  18  is  positive 
definite  for  k,  a  >  0  .  In  particular,  it  is  of  interest  to  select  a  P  which  is 
independent  of  k,  since  x'Px  is  then  a  Liapunov  function  of  the  nonlinear 
system. 

Since  P  is  to  be  a  constant  matrix,  selecting  arbitrarily  p^  =  2 
and  maximizing  Eq.  20,  p^->  =  a  and 


V  =  x' 


a 


(21) 


is  a  CLF  for  the  range  0  <  k  <  k  =  a2  and  indicates  stability  of  the  nonlinear 
system  when  k  =  k  (  x^  )  . 


C.  The  Liapunov  Function  V  =  x'P(k)  x 
The  inability  to  obtain  a  CLF  of  the.  assumed  quadratic  form  (x'Px) 
for  the  linear  problem  in  the  entire  range  k(x^ ) >  0  is  a  serious  drawback 
to  the  method.  The  difficulty  may  arise  due  to  the  restrictive  nature  of  the 
definition  of  a  CLF  ,  i.  e. ,  the  assumption  that  V(x)  is  a  quadratic  form. 

The  off-diagonal  term  ^P^^  ° *  Eq.  17  is  responsible  for  the  upper  bound  on  k. 
Additional  freedom  may  be  obtained  by  permitting  p^  j  to  be  linear  in  k,  i.  e. , 

Pn  =  c  +  vk  (22) 
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constants  . 

Thus,  Q  ( k ) 

of  Eq. 

2  kpl  2 

kp22  ' 

kp22  -  c  - 

vk  +  ap^ 

2ap22 

Since  k  .=  0  lies  in  the  allowable  range,  it  is  necessary  that 

c  =  ap12  . 

With  Eq.  24, 

q(k)  =  det  O  (k) 

becomes 

2 

q(k)  =  -  {kp22  -  vk)  +  4kp12  (ap22-  p12)  . 

If  one  selects 

V  =  P22  » 

then  q  (k)  >  0  implies  no  upper  bound  on  k  as  Eq.  20  did,  since 
q<k)  =  4kp12  (ap22  -  P12  )  . 

Choosing  p22  =  2  and  p^2  =  a 

q (k)  =  4ka2  . 

The  elements  of  the  P  matrix  are 


P11  a  +  2k  >  p12  a  5  p22  ”  2 


The  procedure  followed  so  far  is  to  select  Q  (k)  in  such  a  manner 
as  to  be  least  restrictive  on  the  parameter  k  .  Since  Q  (k)  is  positive 
definite  for  all  k  >  0  ,  P(k)  must  also  be  positive  definite  for  all  k  >  0  ,  since 
it  is  merely  the  solution  of  Eq.  3  for  this  specific  example. 
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It  is  important  to  realize  what  the  modification  of  Eq.  34  has  accomplished 
and  under  what  conditions  such  a  modification  is  possible.  Since  the  linear 
system  is  asymptotically  stable  for  k  >  0  ,  any  Q  (k)  >  0  must  yield  a 
P(k)  >  0  .  In  the  second-order  system  considered  above,  P  (k)  could 
be  written  as  _P^-which  is  independent  of  k—  together  with  a  matrix  which  has 
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all  elements  zero  except  the  element  in  the  first  row  and  first  column  which  is 
2k  .  In  replacing  Eq.  32  by  34  it  is  necessary  that  k  be  a  function  of  x^ 
only  since  the  time  derivative  of  Eqs.  32  and  34  yield  identical  Q's  considered 
as  functions  of  k  .  Further,  for  k  >  0  Eqs.  32  and  34  also  obey 

V(x)  >  x'  Pq  x  .  (36) 

_Pq  is  a  positive-definite  matrix  since  it  is  a  solution  of  Eq.  15  for  k  =  0  . 
Hence,  V  (x)  of  Eq.  34  is  a  positive-definite  function. 

D,  Generalization 

The  results  obtained  in  the  preceding  section  may  be  generalized  to 

the  case  of  a  system  with  several  parameters  kj,  k^ , - k  .  ^or  tke  linear 

system  one  can  attempt  to  find  a  Liapunov  function  of  the  form 

m 

V (x)  =  x'  Pq  x  +  ^  v.  k.x.2  (37) 

i=l 

where  Pq  is  independent  of  the  k^'s  and  w  are  constants  (i=l - m,  m<n). 

For  the  nonlinear  system,  if  k^  =  k^  (x^)  ,  the  Liapunov  function  may  be 
modified  to 

x.  ' 

m  i 

V(x )  =  x1  Pqx  +  2  ^  Vi  k.(u)  u  du  (38) 

i=l  0 

The  total  time  derivative  of  Eqs.  37  and  38  yield  identical  Q's  as  functions 

of  the  k. ' s  . 

- i — 
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m 

V  (x)  >  x'  P  X  +  >  V.  k.*  x.' 

- - -O-  /j  11  1 


(39) 


i=l 


whe  re 


k. 

i 


Max  f  k.  ( x.  )  1 

x.  1  i  i  J 

i 


Min  [  k.  ( x.  )  ] 

x.  L  i  i  J 

l 


v.  <  0 
l 


v.  >  0 

i 


(40) 


Eq.  39  is  clearly  greater  than  zero,  .since  it  is  a  special  case  of  Eq.  37  which 
is  greater  than  zero.  (  j|  x  |j  /  0  ). 

It  is  also  possible  to  make  a  further  generalization  using  LaSalle'  s- 
extension  [2]  of  Liapunov's  theorem.  According  to  LaSalle's  theorem,  it  is 
permissible  for  V  <  0  as  long  as  V  does  not  vanish  identically  on  a 
trajectory.  This  implies  that  Q  may  be  chosen  to  be  only  positive  semi- 


definite.  For  the 

second-order  system  of  Eq.  23,  if 

Q(k)  = 

fo  0 

1 

(41) 

0  4a 

then, 

Pll  8 

c  +  vk  , 

p12  =  °  ,  2a  p22  =  4a 

(42) 

and,  hence, 

c  =  0  ,  ■ 

v  =  ?22’  P22  =  2  • 

(43) 

For  the  nonlinear 

case  we  have  the  Liapunov  function 

*1 

r% 

V (x )  =  4 

a 

\  u  k  (u )  du  +  2x2 

(44) 

and 


V*  (x )  =  -  4a  x^  . 


(45) 
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Hence,  even  for  the  case  when  Q  (  k)  is  only  positive  semi-definite  the 
remarks  of  the  preceding  sections  apply. 

E.  Summary 

The  results  of  this  section  may  be  summarized  as  follows:  Consider 
the  differential  equation  x  =  F  x  which  depends  linearly  on  parameters 

k.  ,  k., - k  .  If  the  solutions  of  the  equations  are  stable  for  k.  <  k.  <3c. 

1.2m  ^  —l—i—i 

1.  V  =  x'  P  (kj  , - km  )  x  is  a  Liapunov  function  if  P 

is  the  solution  of  the  equation 

F'  (k,,---k  );  P  +  P  F  (k. ,  — -k  )=  -  Q(k.,  ---k  )  .  (46) 

2.  If  a  Q  >  0  can.be  selected  such  that  P  is  independent 

of  the  kJs,  V(x)=jx'  Px  is  a  common  Liapunov  function 
in  the  range  k.  <  k^  <  k.  and, consequently,  for  the  nonlinear 
system  where  k.  =  k.  (x)  .  However,  in  most  cases. such  a 
Liapunov  function  can  rarely  be  obtained. 

3.  A  modification  (which  may  be  used  for  the  nonlinear  system) 
is  possible  if  k^.  =  k^  (x^.  ).  In  such  a  case  if  a  Q  >  0  is 

selected  so  that  the  solution  of  Eq.  46  yields  a  P(k^, - k^) 

such  that  p„  is  a  constant  for  i  /  j  and  p„  is  a  linear 
function  of  k^  ,  the  Liapunov  function  for  the  nonlinear 
system  is  as  shown  in  Eq.  38  .  It  is  greater  than  zero  for 

jjx  ||  ^  0  as  it  may  be  bounded  as  in  Eq,  39  by  a  quadratic 

form  coming  from  the  solution  of  Eq.  46  for  a  particular  value 
of  the  constants  k.  in  the  allowable  range. 
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F.  Conclusion 

The  procedure  presented  in  this  section  is  seen  to  be  similar  to  the 

"variable  gradient  method"  [3]  developed  recently.  In  the  conventional 

approach  a  V(x)  >  0  is  selected  as  a  candidate  for  a  Liapunov  function  and 

the  conditions  the  parameters  of  the  system  have  to  satisfy  to  make  V  (x)  <0 

are  determined  to  establish  the  stability  conditions.  In  the  "variable  gradient 

method"  the  gradient  VV  of  V  is  chosen  to  make  V  =  VV'i  negative 
n  (n  —  x  ) 

definite.  — ^  curl  conditions  have  to  be  satisfied  by  VV  to  make  it 
the  gradient  of  some  potential  function  V  (x)  .  One  must  then  insure  that 
V(x)  >  0  for  |!  x  [[  /  0  . 

In  the  approach  presented  in  this  section,  the  nonlinear  problem  is 
related  to  the  corresponding  linear  problem  (when,  all  the  parameters  are 
constant)  for  which  the  stability  conditions  are  known.  By  selecting  a 
Liapunov  function  for  the  linear  problem  and  using  the  procedure  indicated 
previously,  one  is  assured  that  V(x)  is  positive  and  ’O'(x)  negative 
definite  even  for  the  nonlinear  system. 

In  the  following  sections,  the  above  approach  is  applied  to  several 
problems.  Modifications  to  the  procedure  are  also  suggested  -  always 
insuring  that  the  conditions  for  the  stability  of  the  nonlinear  system  reduce 
to  the  Routh-Hurwitz  conditions  when  the  parameters  are  constant. 
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III.  THE  SECOND- ORDER  SYSTEM 
In  this  section^  conditions  are  obtained  for  the  stability  of  certain 
second-order  differential  equations .  The  first  two  results  have  been  discussed 
by  other  authors  [2]  but  are  included  here  for  convenience.  The  conditions 
obtained  for  the  problem  in,  section  C  are,  hopefully  new  and  are  less 
restrictive  than  those  known  to  the  authors. 


A.  x  +  f  ( x ,  x)x  +  g(x)x  =  0 

The  usual  state  vector  representation.for  Eq,  47  with  x  =  x,  is 


2 

-  gx 


1 


fx„ 


(47) 


(48) 


:  For  this  case 


Q  (*,  g)  = 


2gP12 

gp22  -Pu  +  fP12 


gp22  -PH+  fpi2 
2fp22  -  2p12 


(49) 


For-  f  and  g  constant  the  Routh-Hurwitz  conditions  demand  that 
f,  g  >,0  for  stability  or  f,  g  >  0  for  asymptotic  stability.  Hence,  f  =  0  may 
be  assumed  to  be  in  the  allowable  range. 


Q (0, g) 


2gPi2 

Sp22 

g?22  -pll 

-  2P] 

(50) 


c 
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Since  for  g  >  0  Q  (  0  ,  g)  can  never  be  positive  definite,  one  must  take 
Q  (  5  ,  g  )  =  0  so  that 


Pl2  =  0 

Pll  =  p22g  ‘ 

; Since  g  =  g  (x^  )  and  p ^  is  assumed  to  be  a  constant  condition,Eq. 
be  satisfied.  With  Eqs.  51  and  52  ,  Eq.  50  reduces  to 


(51) 

(52) 


52  can 


Q  (f,  g) 


o 


0 


(53) 


so  that  one  may  take  q£2  =  2f  and,  hence,  p^z  =  1 
one  has 


Thus,  for  Eq,  48, 


V  =  2 


A1 

^  ug(u) 


du  +  x. 


(54) 


and 


y  =  -  2f  Uj,  x2).x2 


(55) 


V  may  not  vanish  identically  on  a  trajectory  for  f  >  0  and  g  (x^  )  /  0  for 
Xj  /  0  ,  Thus,  Eq.  54  indicates  asymptotic  stability  in  the  large  for 

(i)  f  (xx  ,  x2)  >  0 

(ii)  g(x1)^'0  for  /  0  (56) 


x 


1 


u  g  (  u )  du 


oo  as 


oo 


(iii)  of  Eq,  56  with  (ii)  demands  that  g(xj  )  >  0  for  x^  /  0  .  [The 
argument  starting  with  Q  (f ,  0  )  was  essentially  given  in  the  previous  section]. 
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B.  x  +  f  (x ) x  +  g  (x ) x  =  0  (57) 

This  is  a  special  case  of  Eq.  47  but  has  been  included  here  to 
indicate  the  importance  of  the  choice  of  state  variables  while  determining 
a  Liapunov  function.  In  a  practical  situation  the  "average"  dissipation  may 
be  more  significant  than  the  instantaneous  dissipation  and  the  Lienard  trans¬ 
formation  [4]  may  be  used.  Defining 


.f 


f  (u)  du  =F(x1  ) 


(58) 


,Eq.  57  reduces  to 


*1  =  *  *  <X1>!X1  +  x2 

x2  =  -  gxl 


(59) 


where 


f  (x1  )  = 


F(xx  ) 

X1 


For  this  case, 


Q(f,.g)  = 

2  (fpn  +  gp12) 

fp12  +  8P22  ’  P11 

A 

fp12  +  gp22  '  P11 

_2p12 

Selecting 

Q  (f,  g)  = 


o 

0 


(60) 


(61) 


(62) 


c- 
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which  is  positive  semi-definite  (or  >  0  )  for  >  0  ,  it  is  found  that  =  0 
and  p^  =  gP22  •  These  conditions  are  compatible  for  p^  =  constant.  Then, 
with  =  Zf'g  ,  one  obtains  -  1  •  Thus, 

X1 

V  =  2  J  ug(u)du  +  X2^  (63) 

0 


and 

V  =  -  2gfxx2  =  -  2gF(xx  )  xl  .  (64) 

Equation  63  indicates  asymptotic  stability  in  the  large  fo* 
a  ^(x  > 

(i)  f(xx  )  =  >  0  for  Xj  4  0 

(ii)  g(x1)>  0  for  Xj  /  0  (65) 

X1 

(iii)  ^  ug(u)  du  -*•  oo  as  [  x^  |  -*•  oo 
0 

With  Eq.  65,  V  will  not  vanish  identically  except  on  the  x  =  0  trajectory. 

Furthermore,  Eq.  65  does  not  demand  that  the  instantaneous  dissipation 
be  >  0  ,  but  only  that  the  "average"  dissipation  be  >  0  .  As  an  example  of 
this,  consider 

,g<*l )  =  1 

f  (xL  )  =  (x1  -  2  )(xl  +  2  )(x1  -/3" )(xl  +/ 3  )  =  xL4  -  7xL2  +  12  . 

Then 


7x, 


3 


+  12 
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which  is  >  0  for  real  .  Notice  f  (x^  )  is  negative  for  -  2<  x^  <  -*JT  and 

yj  <  x^  <  Z  ■ 

C.  Second-Order  Nonlinear  Feedback  System 


J - .  V_ 

J — - — 

— — 1 

r  — 

- 

f 

g 

h 

_ , 

Figure  3,  Second-Order  Nonlinear  Feedback  System 


The  system  of  Fig. 

x 

* 

X 


3  may  be  described  by  the  set  of  equations 

1  =  +  x2 
2  =  ‘h(xl)xl  -  f<x2)x2  ‘ 


(66) 


This  may  be  reduced  to  the  single  second-order  differential  equation 


x  + 


g  (  X  )  +  g'  (  X  )  X  +  ^1  +  -g  ^ x  ~X—  j  f  [x  +  g  (x )  xj 


x  +  h(x)x=0  (67) 


where  x  =  x^  . 

Eq.  67  is  the  standard  form  of  Eq.  47.  There  are  two  immediate 
objections  to  applying  the  results  of  A  .  A  perhaps  minor  objection  is  that 
Eq.  67  contains  a  derivative  of  the  function  g  .  It  would  be  desirable  not 
to  demand  that  the  nonlinearities  be  differentiable  functions.  A  more  serious 
objection  and  one  in  fact  that  prohibits  the  application  of  A  ,  is  that  the 
"dissipation"  of  Eq.  67  is  not  of  constant  sign  in  the  vicinity  of  x  =  0  . 
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Examining  Eq.  67.  for  the  case  of  f,  g,  h  constant,  one  obtains 

ii  +  (f  +  g)x  +  (fg+h)x=0  .  (68) 

The  system  is  thus  stable  for 

( i )  fg  +  h  >  0 

(69) 

(ii)  f  +  g  >  0 

As  has  been  previously  mentioned,  one  wishes  to  select  the  Q(f,  g,  h)  for 
the  nonlinear  system  which,  will  be  no  more  restrictive  than  Eq.  69  when 
applied  to  a  linear  system. 

Q(f,  g,  h)  = 


|2(png  +p12h  p22h +pl2(f+g)  -  pn 

p22h  +  p12^  +  &)  '  P11  2(p22f  "  p12  ^ 

As  in  Eq.  23,  the  off-diagonal  elements  of  Eq.  70  will  introduce  undesirable 
restrictions  on  f,  g,  and  h  .  For  f,  g,  h  constant,  one  could  select 


(70) 


Q(f,g,h) 


*11 

0 


(71) 


and,  hence,  obtain 


p12  =  P22 f 

pll  =  p22h  +pi2(f+g) 

*11  =  2  <Pl  1  e  +p12h) 


(72) 


If 


Ml 


4(f  +  g)  (fg  +h) 


(73) 
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then  the  conditions  of  Eq.  69  are  clearly  displayed.  WithEqs.  72  ahd  73  , 
one  finds 

pu  =  2h  +  2f  (  f  +  g  ) 

P12  =  2f  (74) 

?22  =  2 


For  the  nonlinear  problem,  however,  since  f  =  f  (x^),  one  is  unable 
to  select  pji  and  p^  as  Eq.  74  •  Realizing  that  Eq.  69  implies  that 
for  g>  0  any.  f  greater  than  some  minimum  is  satisfactory,  one  may  try 


pn  =  2h  +  2f_  (_f  +  g) 

Pl2  =  2L 

PZ2  =  2 

where 

£  =  Min  [  f  ( Xg  )  ■]  . 

x2 

Notice  that  even  for  f,  g,h  constants  and  g  <  0  and  h  >  0  , 
Eq.  69  is  not  permitted  to  increase  without  bound. 

With  Eq.  75  ,  70  becomes 


Q 


2  (h  +ig)  (f  +  g  ) 
f  (f  -  f) 


whicH  i®  at  least  >  0  for 


f  (f  -  f) 

2  (f  -  f) 


(h  +  fg).(f  +  g) 


.  iz  (£  -  J)  . 
>  — ? -  £ 


0 


f  from 


(7  5) 


(76) 


(77) 


(78) 
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since  f  >  £  .  If  one  defines 


f  =  Max  [f  (x2)] 
x2 


then  the  system  is  asymptotically  stable  in  the  large  for 
(i)  f  g  (x1  )  +  h  (x1  )  >  0 

2 

(ii)  tf+  g  (x^]  [fg  (xx).+  h  (xj)  ]  >  >  0 


and 


and 


l 

(iii)  ^  [f_g(u).  +  h  (  u  )]udy.  -+  oo  as  |xj  (  <x> 

*0 

X1 

V  =  2  ^  [  f  g  (  u)  +  h  ( u)]  u_du  +  {fxl  +  x-,)2  . 


(79) 


(80) 


(81) 


Eq.  80  reduces  to  the  conditions  of  Eq.  69  for  f  constant  even  with 
g  =  g  (xi  )  and  h  =  h  (x^  )  and  also  permits  a  negative  g  ( x^  )  .  (ii)  of 
Eq.  80  may  be  used  to  determine  the  maximum  value  f  ( x2  )  may  assume 
while  still  being  guaranteed  stability  by  Eq.  81  .  It  is  clear  that  the 
conditions  of  Eq,  69  are  necessary  to  allow  for  the  validity  of  Eq.  80  . 
Further,  notice  that  it  was  not  necessary  to  introduce  the  derivatives  of  any 
of  the  functions  through  a  Jacobian  matrix  [5].  This  was  partially  due  to 
the  form  of  F  ,  i.  e.  ,  F. .  =  F. .  (x. )  only. 

As  an  example  of  the  above,  take 

-1*^  | 

f  (xl)  =  2  +  ae  ,  a>0 

2 

"X1 

g  (xj  )•  =  -  e 

and 

h  (xx  ).  =  6  +  2  |  |x1 1 
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Then  f  =  2  and  f  =  2  +•  a  for  a  >  0  .  It  is  desired  to  find  how  large  a 
may  be  taken.  Now  { j  g  +  h )  >  (-2  +  6)=4>0  and  _f  +  g>2-l  =  l>0  . 
Thus, 


(f  g  +  h  )  (f  +  g  )  >  4 

One  may  select,  therefore,  from  Eq.  80 


jjjf  -1) 

4 

and,  hence, 

f  =  _f  +3 

so  that  stability  is  assured  for  0  <  a  <  3 


3  , 

5  , 


IV.  THIRD-ORDER  SYSTEM 

In  this  section,  conditions  for  the  stability  of  the  differential  equation 
x  +  f(x,x,x)x+-g(x,x)*+h(x)x  =  0  (82) 

are  considered.  This  equation  can  be  represented  in  the  standard  state  vector 
form  as 

*1  =  x2 

*2  =  x3  (83> 

*3  =  "  h(x1)x1  -  g(xj,  x2)x2  -  f(xj,  x2,  x3)x3 
1  ' 


where  x  =  x 
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When  f,  g  and  h  are  constants,  the  Routh-Hurwitz  conditions  for 
stability  are 

(i)  f ,  g,  h  >  0 

(84) 

(ii)  f  g  -  h  >0  . 

The  conditions  for  stability  of  the  nonlinear  problems  to  be  derived  below 
are  such  that  they  reduce  to  Eq.  84  when  f ,  g  and  h  are  constant. 


A.  x  +  f  ( x  )  x  +  g(  x )  x  +  x  =  0 


Using  the  above  notation, 
Q  (f,  g,  h.)  =  2P13h 


P13g+P23h-Pll 


p33h  +  p13f'p12 


p13«  +  p23h'pir  2(p23S'p12)  p33*+p23f-p22"p13 

p33h  +  p13f‘p12  p33g+p23f-p22'p13  2^p33f'p23) 


The  problem  at  this  stage  is  to  determine  Q  so  that  it  contains  in  a 


readily  recognizable  form  the  conditions  of  Eq.  84  if  f,  g,  h  were  constants. 
One  may  thus  proceed  as  in  Section  III-C  and  attempt  to  employ  a  diagonal 
Q.  Noticing  that  q23  =  q32  contains  f  ( x2  ) ,  g(x2),  and  p22  ,  it  is 
possible  to  make  these  terms  zero.  Thus,  one  tries 
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Solving  Eq.  85  with  86  , 

Pll  =  ?23h 
Pl2  =  P33 h 
Pl3  =  ° 

P22  =  P33S  +  P23f 

q22  =  2(P23«  -  P33h) 

q33  =  2  (p33f  -  p23  ) 

.Selecting 

q22  =  2  (mg  -  h) 

and 

q33  =  2  (f  -  m) 


where  m  is  a  constant, 


and 


P23  = 
P33 


then 

m 

1 


Since  Q  must  be  positive  definite,  it  is  necessary  that 
(i)  m  >  0 

and 

(ii)  m  g  -  h  >  0 


(87) 


(88) 


(89) 


(90) 


where 
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[Equality  is  not  permitted  for  both  (i)  and  (ii)  of  Eq,  90]  .  m  may  be  selected 
to  maximize 

q  .=  (f_  -  m)  (mg  -  h)  (93) 

or 

Ig  +  h 

(94) 


m  = 


2& 


With  Eq.  94  ,  Eq.  90  reduces  to 

£g-  h 


(i) 


2g 


>  0 


Ig.-  h 
(ii)  - i —  >  0 

Thus,  for  Q  >  0  ,  it  is  necessary  that 
(i)  g  >  0 
(ii)  ,lj;  -  h  >  0 

For  (ii)  of  Eq.  96  to  hold  for  h  =  0  ,  one  must  also  have  I  >  0  .  Taking 
h  >  0  ,  it  is  assured  that  m  >  0  which  is  required  by  =  mh  . 

In  summary,  the  system  is  asymptotically  stable  in  the  large  for 

(i)  1,1,  h  >  0 

(ii)  I  £  -  h  >  0 

where  x^ 

V  =  mhxj2  +  2hXj  x.,  +  2^  [g(u)+mf(  u)]  u  du 


(95) 


(96) 


(97) 


+  2mXj,  x^  +  x^* 


(98) 
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and  m  is  given  by,  Eq.  94  .  Further, 


where 


* 

P  = 


V  >  x'  P  x 

mh  h  0 

h  ( j[  +  mf)  m 

Oml 


(99) 


(100) 


which  is  positive  definite.  With  Eq.  97,  V  may  not  vanish  identically  ex¬ 
cept  on  the  x  =  0  trajectory.  For  f  and  g  constant,  conditions  Eq.  97  are 
identical  to  those  of  Eq.  84. 


B.  x'+  f(x)x  +  g(x)  x  +  h(x)  x  =  0 


f 


This  case  is  to  be  identical  to  that  of  A  ,  except  here  h  is  not  a 

constant  but  h  =  h  (x^  )  .  If  Eq.  98  were  independent  of  h,  then  it  would  also 

be  a  Liapunov  function  for  this -case.  However,  Eq.  98  is  a,  function  of  h 

2 

but  fortunately  only  in  what  corresponds  to  the  P^xi  term  and  the 

2p i2xix2  term  (taking  m  as  constant).  A  modification  of  these  terms  may 

be  attempted  in  a  manner  that,  will  not  drastically  alter  the  Q  of  A  . 

By  examining  Eqs.  8  5  and  87  ,  it  is  clear  that  the  proper  generalization 
2 

of  the  term  is  to  replace  it  as  before  by 


2 


I 


Pu  ( u)  u  du 


0 

since  it  came  from  p^  =  ?23^  ’  ^or  t*1*8  ca8e  then,  one  might  try  this  as 
the  only,  change  taking 
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x2 


V  =  2m  (  u )  u  du  +  2h  ( Xj  )  Xj  x2  +  2  |*  [g(u)  +  m£(u)  ]  u  du 


(101) 


+  2m  x2  +  x^ 


It  is  clear  that  m  has  to  be  redefined  from  Eq.  94  for  it  is  required  to 
be  constant. 

Trying  Eq,  101  ,  one  obtains 


where 

Q  = 


0 

0 

0 


V  =  -x'Qx 

0 
0 

2  (  f  -  m  ) 


0 

2  (mg  -  H') 
0 


(102) 


(103) 


and 


H(x1  )  =  h(xx  )  Xj  , 


H' 


dH 

dxj 


(104) 

(105) 


Thus,  for  Q  >  0  even  in  the  most  restrictive  case,  it  is  necessary  that 
(i)  X  -  m  >  0 

and  _  (106) 

(ii)  m_£  -  "TT"  >  0 

where 


7F  =  Max  [H'  (xx  )  ]  .  (107) 

Here  again  equality  is  not  permitted  for  both  (i)  and  (ii)  of  Eq-106.  If  m  is 
selected  as  in  A  to  maximize  the  product  of  (i)  and  (ii)  of  Eq.  106,  then 


+  W 

m  =  .  - 

2£ 


(108) 
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With  Eqs.  106  and  108 
(i)  g_  >  0 
(ii)  f  g  -  H1  >  0  . 


following  A  ,  for  Q  >  0  ,  it  is  necessary  that 

(109) 


If  (ii)  of  Eq.  109  is  to-hold  for  h  =  0  ,  one  must  have  £>  0  .  Here  H'  >  0 
guarantees  an  m  >  0  . 

With  m  defined  by  Eq.  106  and  taking  V  of  Eq.  101  as  the  Liapunov 
function,  Xj 

V  >  2m  y  H(u),du  +  2H(x^  )x2  +  gx^2  +  mfx22  +  2mx^  x3  +  x32  ( I  ip) 

0 

or 


V  > 


H(x.) 

- - -  + 

.vT 


+ 


[x2»x3] 


mf 

— 1 
m 

- 

X-. 

m 

i 

2 

_x3 

J. 

+  ~  y  [m  g-H'(u)]h(u)  u  du  • 

~  0 

(111) 


■With  Eq.  109  ,  m  >  0  ,  and  h  (x^  >  0  for  4  0  ,f'V.  of  Eq.  Ill  is  >  0 
for  ||  x  ||  4  0  . 


(iv) 

only 


To  summarize,  the  system  is 
(i)  h  (x^  )  >  0  for  Xj  /  0 
(ii)  £,  g ,  IT  >  0 

(iii)  fg  -  H7"  >  0 
X1 

(iv)  jj*  H(  u )  du  oo  as  |  x^  j 
0 


asymptotically  stable  in  the  large  for 


(112) 


CO 


of  Eq.  112  insures  that  V-*oo  as  ||x||  -*■  oo  .  Also  notice  that  V  may 
vanish  identically  on  the  x  =  0  trajectory. 
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C.  The  General  Case  x  +  f  (x,  x ,  x )  x  +  g(  x,  x )  x  +  h  (x )  x  =  0  , 

It  is  desirable  that  the  results  of  this  section  include  as  a  special 
case  the  results  of  the  sections  A  and  B  .  This  implies  that  a  Liapunov 
function  must  be  obtained  which  would  reduce  identically  to  that  of  Eq,  101 
for  f  =  f  (x^)  and  g  =  g  (x^)  .  To  this  end,  one  anticipates  a  redefinition  of 
m  of  Eq.  108  and  a  change  in  the  integral  involving  f  and  g  .  Thus  try 
X1  x2 

V  =  2m  ^  H  (u)du  +  2H  (xj  )  +  2  J  [  g  (u)  +  mf  (u)  ]  u  du 


0 


0 


+  2m  x^  x^ 


+  x_ 


(113) 


where  g  ( x^ )  and  f  (x^)  are  to  be  determined  in  a  manner  to  reduce 
identically  to  g  and  f  when  they  are  functions  of  x.,  only. 

Now 

ir  =  -  x'Qx 

where 

0  0 
2  (mg  -  H'  )  (g  -  g)  +  m  (f  -  $ 

(g-g)  +  m  (f  -  f')  2  (f  -  m) 

To  insure  that  Q  >  0  ,  one  demands  that 

(i)  mg  -  H>  >  0 

(ii)  f  -  m  >  0 

(iii)  q  =  4  (mg  -  H1 )  (f  -  m)  -  Max  [(g-g)+m(f-f)]2>0 


(114) 


(115) 


and 


(116) 
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where  here  the  =  means  the  absolute  minimum  of  the  particular  function, and 
H'  was  defined  in  Eq.107.  If  Eq.116'  is  satisfied,  then  Q  >  0  . 

(iii)  of  Eq.  116  is  the  requirement  that  was  not  present  in  the  previous 
sections.  This  condition  will  be  the  least  restrictive  if  one  selects  g  and  f 
so  as  to  minimize  Max  [(g  -  £)  +  m  (f  -  f)  ]2  .  With  the  stipulation  on  g  and 

A 


f  of  the  first  paragraph  of  C  ,  it  is  natural  to  make  the  following  definitions.: 


or 


f  ( x?  )  =  Max  [f  ( x,  ,  x 

2  >  x3  )  ] 

xp  x3 

f  (x2 )  =  Min[f(x^,x 

2  >  x3 )  1 

xl>x3 

A  f  (*?)  +  UX?) 

f (x2)  = 

2 

f  =  Max 

Max 

f(x  ,  X 

x2 

X1  >  x3 

A  i 

T(x2  ) 

-  f(x2r 

f  =  Max 

7 

x2 

Ld 

_ 

(117) 


-  Min  [i  (  x2  )  ] 
X2 


The  same  notation  is  adopted  for  g  (x-^  ,  x^  )  ,  e.  g.  g  (  x-> )  =  Max  [g(x, ,  x,)], 

^  _  Xj  * 

to  define  g  (x2)  ,  g  (x2)  ,  g  (x2  )  ,  g,  and  g;  .  H'  has  been  defined  previously. 
With  the  above,  (iii)  of  Eq,  116  requires  that 
q  =  4  (m£  -  H'  )  (X.  -  m)  -  ( g"  +  mT  )2  >  0 


(118) 
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If  now  m  is  selected  to  maximize  q  of  Eq.  118  ,  it  is  found  that 


m 


2  lg  +  2H' 
f  +  4g 


-  f  g 


With  m  taken  as  in  Eq.  119  ,  it  is  found  that  q  >  0  requires 


fg  -  H'  >  € 

where 

€  =  \/  T  ^  £  H'  +  T'  g'  (£g  +  H<  )  +  g'  Z  g 

To  summarize, the  following  conditions  insure  that  Eq.  82  is 
asymptotically  stable  in  the  large: 

(i)  h  ( )  >  0  for  Xj  /  0 

(ii) i,  g,  H1  >  0 

(iii)  -  H'  >  e  >  0 
X1 

(iv)  f  H  (  u)  du  -*■  oo  as  |x^  |  -*■  oo 


(119) 


(120) 


(121) 


(122) 


where  e  is  defined  in  Eq,  120  and  the  Liapunov  function  is  defined  in 
Eq.  113  with  m  given  by  Eq.  119  .  With  Eq.  122  ,  one  can  show  as  in  B 
that  V  >  0  for  j  |xj  |  /  0  since  Eq.  122  implies  m  >  0  .  Further,  V  may 
only  vanish  identically  on  the  x  =  0  trajectory. 

The  conditions  for  stability  of  Eq.  122  and  e  in  Eq.  121  may  appear 
complicated.  The  following  four  cases  ,  however,  indicate  that  these 
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conditions  may  be  reduced  to  simpler  forms  for  . specific  types  of  npnlinearities 
and  do  contain  the  results  of  A  and  B  . 

1.  Take  f,  g,  and  h  as  constants.  Then  f  =  g  =  e  =  0  ,  and-Eq,  122 
reduces  to  the  ilouth-Hurwitz  conditions  of  Eq.  84. 

2.  Take  f  constant,  g  =  g  (x),  and  h  =  h (x)  .  Then  f  =  g  =  e  *  0  , 
and  Eq,  122  becomes  Barbashin's  result  [6]  . 

3.  Take  f  =  f  (x)  ,  g  =  g  (x)  ,  and  h  =  h  (x)  .  Then  e  =  0  and  Eq.  122 
is  just  Eq,  112  of  B  .  For  h  =  constant,  Eq.  122  becomes  Eq,  9 7 
of  A  . 

4.  Cases  2  and  3  considered  above  require  conditions  similar  to 
those  of  a  linear  system.  Ifg  =  g(x),f  =  f(x,.x),  and 

h  =  h  (x)  an  e  may  be  determined  satisfying  (iii)  of  Eq.  122. 
Since,  in  this  case,  g”  =  Oandji  =  f, 


z7=  r-i=  — Si 

■  yw 

and  the  condition  reduces  to 


£<£  (x ,  x)  <T=  £ 


+ 


2e 

yw 


[Numerical  example:  assume  ^  =  I  ,  g  =  7  ,  and  "IF  =  4  so  that 
Jg  =  IF  =  3  >  0  .  Taking  e  =  2.  99  it  is  seen  that  1  <f(x,'x)<  3.  99 
insures  stability.  ] 


I 
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Consider  the  .following  example  which  indicates  an  application  of 
the  most  general  results.  Assume 

-[x.^  +  x3  +  x3] 

f  (x^  ,  )  =  ae  .'  +  1  +  a 


x2  +  1 


g  (xx  ,  x2)  =  6  +  be 


r  6  ,  2n 

■[*l  +  x2  ^ 


2,b 


x2  +  2 


"X 

h(xL  )  =  1  -  [1  -  (1  +  xx  )  e  1  ] 


for  Xj  >  0  and  h  (  -  Xj)  =  h  (x^) 

It  is  desired  to  determine  tHe  range  of  a  and  b  so  that  stability  is  assured. 
Using  the  previous  definitions 
-x2 


f(x2)=ae  *  +  1  + 
f  (x2)  =  1 


f  (x2)  =-r 


f  =  Max  -g- 
x2 


2  ,  1 

e  + 


x2  |  +  1 


+  1 


2  , 
e  + 


x2 


+  1 


=  a 


f  =  1 


g  (x2)  =  6  +  be  2  +  — 


x2  +  2 
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Also 


and 


and 


g(x2)  =  6 


g  (-3C,)  =6  +-J- 


'X2  , 
e  + 


2 

x-,^  +  2 


Max  -j- 
x2 


2 

e  + 


“T~“ 

x2  +  2 


£ 


=  6 


Hlxjlsxj  -  3  [1  -  (1  +  xx)  e  1  ] 


H(  -  Xj)=  -  H(xj) 
H'  (Xj )  =  1  -  3  Xj  e 


H'  (  -  x2  )  =  H'  (Xj  ) 

EF  =  1  „ 


for  x}  >  0 


for  Xj  >  0 


Thus,  _f  £  -  TF  =  5  >  0  and  one  may  select  e  <  5  .  Then  with  Eq.  121  , 
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The  region  of  allowable  a  and  b  is  shown  in  Fig.  4. 


u 


Figure  4.  Region  of  Permissible  Values  of  a  and  b 


For  example,  one  may  take 


1. 

a  =  0 

and  b  < - 

S* 

2. 

o 

u 

rO 

and  a  <  5 

3. 

a  =  b 

<  5 

/I4 

All  the  conditions  of  Eq.  122  have  been  obeyed.  Notice  that  H'(xj  )  is 
negative  for  certain  values  x^  ,  e.  g. ,  |  Xj  f  =  1  . 
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D.  x  +  fx  +  g(x)x  +  hx  =  0  .  (123) 

In  this  section,  f  and  g  are  constants.  Of  course,  Eq.  123 
may  be  considered  as  a  special  case  of  C  .  However,  less  stringent 
conditions  will  be  derived  for  stability  than  those  given  by  Eq.  122  for  this 
special  case. 

The  following  state  vector  form  is  taken  which  may  be  considered  as 
a  logical  extension  to  Eq.  123  of  the  Lienard  transformation  of  Section  III  -  B 
for  a  second-order  system: 

Yi  =  y3 

y2  =  -  hy1  (124) 

y3  =  -  G(yj  )■  +  y2  -  fy3 

where 

yx  =  x 

and  y^ 

Gfy^  =  J  g(u)da  .  (125) 

0 


Q 


Proceeding  as  in  A  , 


2[hpi2+-^-  pi33  hP22+  ■~P23-P13 

Q 

hp22  +  y^p23"p13  '2p23 


hp23+  ^p33+fp13‘pll 
fp23"pl  2  ”  p33 


hp23+yi  p33+fp13"pll  fp23"p12"p33 


2  (fP33_Pl3) 


(126) 
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Notice  that  for  g  constant, 


G(y1 ) 

?1 


One  may  select 

Q  = 


so  that 


2(f  —  -  h) 
*1 

0 

0 


_  ,2  ,  G 
PH  -  f  + 


?1 


Pl2 


=  -  1 


p22  =  h 
p33  *  1 


p13  =  f 
p23  =  0 


and,  hence, 

V  =  f2y12  +  2j  G(u)du  -  2yi  y2  +  ^-y22  +  Zfy^  +  Y $ 


*1 

0 


or 


(127) 


(128) 


(129) 


l 

=  t  I[f  ^r1  -  h]  udu  +  4  {yz  -  t~  y\)Z  +  (iy\  +  ^)2  • 


(130) 


Thus,  to  summarize,  Eq.  123  is  asymptotically  stable  in  the  large  for 

°(yi)  . 

(i)  __ -  ,  h  >  0  for  y,  /  0 

y  1 

<%!> 

(ii)  f  - —  -  h  >  0  (131) 

y\ 

?i 

(iii)  ^  [  f  -  h  ]  u  du  -*•  oo  as  |  y^  j  .-*>  ao 
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Further,  V  >  0  .for  |  |jrJ  |  y  0  and  that  with  (iii)  of  Eq.  131,  V -►  oo  as 
i  III  I  "*■  cd  .  Also,  with  Eq.  131  ,  V  may  only  vanish  identically  on  the 
y  =  0  trajectory. 

G(yi) 

It  is  required  by  Eq.  131  that  only  — — -  be  >  0  and  not  that 

g(yL  );>  0  .  As  an  example  of  this,  consider 
f  =  36 


g(yx)  =  (yx  -  2)  (y2  +  2)  (y2  -J%)  (yj  +yi)  =  -  7yx2  +  12 

h  =  186  . 


Then 


(yi} 

yi 


7y, 


+  12 


G  ( y  i  )  i  o  n 

which  is  >  0  for  y^  real.  Further,  the  minimum  of  — —  is  -jtj- 


Thus,  f- 


tYl> 

y\ 


187 


h  >  f  —jg—  -  h  =  1  >  0  .  All  the  conditions  of  Eq.  131 


are,  hence,  satisfied. 


TR403 


-39- 


V.  CONCLUSIONS 

The  various  results  presented  throughout  this  report  were  derived 
by  considering  a  CLF  (or  a  modification  of  the  CL-F  )  for  a  related  linear 
system.  All  the  results  for  the  nonlinear  equations  reduce  to  the  usual 
Routh-Hurwitz  conditions  when  the  nonlinearities  are  assumed  constant. 
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